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1. INTRODUCTION 



The precision determination of the Cabibbo-Kobayashi-Maskawa (CKM) matrix element 
\Vah\ is of particular importance to test the Yukawa sector of the Standard Model (SM). 
More precisely, it's needed to test the mechanism which explains the occurrence of Charge 
and Parity (CP) violating effects in weak decays. Such effects are caused by the presence of 
an irreducible complex phase in the unitary 3x3 CKM matrix which can be illustrated by 
using one of the triangle equation derived from the unitarity constraint of the CKM matrix, 

1. e. Kid Kib + ^cd Kb + ^td Kb — 0- Given the proper normalization of the sides of the triangle 
the complex phase corresponds to the apex of the so-called unitary triangle in the complex 
plane, what is illustrated in Fig. [1} The left side of the unitary triangle is proportional 
to the absolute value of V^b, which can be combined with the direct measurements of the 
angles of the triangle and the other CKM matrix elements to test the predicted unitarity of 
the SM CKM matrix. Such tests were performed e.g. by the authors of Refs. [52] and 
and found an excellent agreement with unitarity within todays experimental precision. 

The CKM matrix element |Vub| can be determined from a multitude of weak decays 
which involve either inclusive or exclusive final states and exhibit different experimental or 
theoretical challenges. Fig. 121 illustrates the leptonic, hadronic, and semileptonic transitions 
whose partial or total decay rates are proportional to ^ |Vub| ■ In this presentation I will 
focus on the determination of |Kib| from semileptonic decays, which have some definite ad- 
vantages over the determination of |Kib| from leptonic, e.g. through B — )■ rz/,-, and hadronic 
decays: The leptonic determination is experimentally challenging, and the prediction of the 
hadronic decay rate involves complicated non-perturbative and perturbative corrections due 
to the strong interaction processes between the two hadronic final states. The study of 
semileptonic decays offers some middle ground between experimental and theoretical chal- 
lenges: the presence of a high-energetic lepton offers a good discriminator from other weak 
decays, and the determination of the decay rate is simplified due to the factorization of the 
matrix element into a hadronic and leptonic current. In practice, however, the presence of 
the much more abundant semileptonic 6 — >■ c transition complicates things and one neither 
has to identify one particular exclusive final state (e.g. one pion), restrict the measurement 
to regions of phase-space where charmed semileptonic decays are kinematically suppressed 
or forbidden, or use data mining algorithms building on multivariate methods which can 
separate charmed and charmless semileptonic transitions. The presented results in this talk 
determine |Vub| either via the reconstruction of an exclusive final state, or measure the fully 
inclusive (partial) decay rate. 

2. INCLUSIVE DETERMINATION OF \Vub\ 

For inclusive B — )■ Xu I i^i decays, where Xu denotes a hadronic system consisting of one 
or more mesons with at least one u quark produced in the weak h ^ u quark transition, 
a prediction for the total decay rate can be readily obtained in the heavy quark expansion 
(HQE) of the Standard Model Lagrangian. The semileptonic transition can be expressed 
as a sum of local operators using an operator product expansion (OPE), cf. e.g. Ref. [37] 
for a review. Unfortunately the dominant B — )■ Xf-lvi decays do not allow the direct 
measurement of the total decay rate. Many analyses are thus forced thus to measure the 
partial branching fraction in a region of phase-space, where the background from 6 — )• c 
transitions is considerably suppressed or kinematically forbidden. Although in these regions 
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Fig. 1 The unitary triangle for V^d V^b + ^cd + = is shown. The illustration was 

taken from Ref. 1461. 




Fig. 2 The Feynman graphs for the leptonic, semileptonic, and hadronic decays of a i? meson 

are shown. The illustrations were taken from Ref. |46j . 



a clear separation from charmed semileptonic decays emerge, the OPE breaks down and non- 
perturbative and perturbative correction to the partial decay rate become very important. 
In particular, the Fermi motion of the h quark inside the B meson becomes a relevant factor 
which e.g. determines the shape of the lepton momentum spectrum near the endpoint or 
the hadronic invariant mass distribution, m^, for small values of mx- This behaviour of 
the OPE can be remedied by summing the most divergent contributions of the expansion 
into a single function, the so-called shape function. In practice, the functional form of the 
shape function is unknown. Its moments, however, can be related to the moments of the 
measured mass and lepton spectra of 6 — )■ c decays and observables from 6 — )■ 57 decays, 
cf. e.g. Refs. [38j. Thus the predictions of the partial decay rates of Refs. [H [SI |35l HI] 
employ a variety of model functions to describe the shape function. The first few moments of 
these model function are matched to the experimental available input, and are used under 
the premise that the precise functional form of the shape function only has a negligible 
impact on the resulting partial decay rate. Experimentally, the use of multivariate methods 
allowed reducing the dependence of the shape function by measuring a large fraction of the 
total phase space. In this talk the recent measurements of Refs. [S^ is reviewed whose 
authors successfully employ such a multivariate technique to measure about ^ 90% of the 
B — )■ Xu I vi phase-spac. Furthermore, I review the recent analysis of Ref. [42j, whose authors 
use a veto based approach and employ phase-space cuts to reject charmed background. 

A. Study of B ^ X^lPi from Belle: Ref. [53] 

The authors of Ref. [53] measure the partial B — t- X„ I vi branching fraction with a non- 
precedented coverage of the allowed B — )■ Xu I vi phase-space. In total 657 x 10^ BB pairs. 
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measured at the KEK-II B factory with the Belle detector, were analyzed in a tagged ap- 
proach, i.e. one of the decaying B mesons is completely reconstructed to infer the kinematics 
and flavor for the recoiling B meson. This so-called tag or reco i?-meson is reconstructed 
using a multitude of exclusive hadronic decay modes. For each such tag B meson candidate 
the beam-energy substituted mass 

mbc/ES = \/s/4:- \p^f , (1) 

and the energy difference 

AE = E%- (2) 

is reconstructed, where {E^,p^) is the four-momentum of the tagged 5-meson candidate 
in the T(45') rest frame, and ^/s denotes the beam energy. For correctly reconstructed 
reco S- meson candidates one expects rrihc/ES to peak at the i?-meson mass, and AE at 
zero. The primary vertex of the i?-meson reco candidates is determined using a vertex fit 
and in if more than one candidate is present in an event, only the one with the highest 
significance with respect to the quality of the fit is retained. Further, a quality cut on the 
value of the vertex fit is implemented. The authors require a reco candidate 5-meson 
to lie within 5.27 GeV < nibc/ES < 5.29 GeV and -0.05 GeV < AE < 0.05 GeV. The 
shape of combinatorial background in the given rrihc/ES and AE window is estimated from 
Monte Carlo (MC) simulations and subtracted from the candidates. Contributions from 
non-BB continuum background is subtracted by a off-resonance data sample which was 
taken below the T(45') resonance scaled by the integrated on- and off- resonance luminosity 
ratio. Electron and muon candidates from the recoiling other B meson are required to pass 
angular acceptance cuts and to be identified by a particle identification algorithm. Due to 
the reconstructed 5-meson four-momentum from the tag side, the lepton three momentum 
can be boosted into the rest frame of the recoiling i?-meson. The hadronic Xu system 
associated with the B — t- Xu I i^i decay is reconstructed from charged tracks and energy 
depositions in the calorimeter that are not associated with the tagged side of the event. The 
further B — )■ X„ / i?i signal selection is based on a non-linear multivariate boosted decision 
tree (BDT), cf. Ref. |^38j. The BDT incorporates a total of 17 discriminating variables to 
create a single classifier, which separate the B — )■ X^ I vi semileptonic decays from other 
background. These discriminative values include non-exhaustively kinematic quantities, the 
number of identified kaons in an event, rribc/ES, the absolute value of the net charge of the 
event, and the track multiplicity, which tends to be higher for recoil candidates originating 
from 6 — )■ c transitions. Further, if a low-momentum pion is present in a given event it is 
associated with a potential strong D* — ^ Dn decay. Due to the low momentum transfer of 
the D* to the D system, the pion three-momentum carries almost all of the three-momentum 
of the D* and can be used to infer the D* four-momentum. This can be used to calculate 
the missing mass squared associated with a. B ^ D* I ui decay, i.e. 

^LssD* = (Pb' - Vd' -Pif , (3) 

where pb' is the four- momentum of the recoiling signal 5-meson inferred from the tagged 
side, Pd* denotes the four-momentum inferred from the slow pion, and pi is the four- 
momentum of the lepton candidate. In case of a true B ^ D* I vi decay, Eq. [3] is expected 
to peak at zero. Using these discriminators as input for the BDT a good signal to back- 
ground separation can be obtained down to a lower lepton momentum cut in the 5-meson 
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Fig. 3 



The one-dimensional projections of the fitted mx and q'^ distributions of 

Ref. [53] are shown. 



rest frame, of 1 GeV/c^. This allows to measure the partial branching fraction of of 

about 90% of the allowed phase-space of the B — )■ I i^i transition. 

The candidates passing the preselection and the selection of the BDT classifier are ana- 
lyzed in a two-dimensional fit in the reconstructed hadronic invariant mass of the X„ system, 
mx, and the four- momentum transfer squared of the fi- meson to the system, i.e. in the 
T(4iS') rest frame 

q'={iV~s,0)-p*B-Px)\ (4) 

where p*^ denotes the four-momentum of the tagged B meson, and p*x is the four-momentum 
of the reconstructed hadronic Xu system. The free parameters in the two-dimensional fit are 
the yields of B ^ Xu luu B ^ Xclvu and other backgrounds caused mainly by secondary 
and misidentified leptons. 

The B — )■ Xu I vi contributions are modelled using a hybrid mix of inclusive and exclusive 
contributions. Resonant contributions from B ^ n I ui, 5— )■ plui, B ^ ulvi are modelled 
using the form factor predictions of Ref. [HI [T3] . The resonant decays into B rjlui and 
B ^ 1]' lui are modelled using the quark-potential model of Ref. . The corresponding ex- 
clusive branching fractions are fixed at the measured world-averages calculated by Ref. [TB] . 
The non-resonant contributions are modelled using the shape function parametrization of 
Ref. [33] . The hybrid MC is required to match the moments of the and mx distributions 
predicted by Ref. [31] . 

The B — )■ Xclvi contributions are modelled as the sum of the exclusive B — )■ Dli?i, 
B ^ D* Ivi, B ^ D** I Pi, and non-resonant B — )■ D^*^ irlui decays modelled after Refs. [25], 
|36l Us]. The non-perturbative shape parameters of the B — )• D^*^ I ui decays are fixed at the 
measured world average of Ref. [TB] , and the branching fractions are fixed at the measured 
values of Ref. [3]. The narrow B — D** / z/; and non-resonant B — t- D^*^nli'i branching 
fractions are fixed at the averaged values of Ref. [TB] . The broad B — )■ D** I ui branching 
fractions are fixed to fill the gap between the inclusive and summed exclusive branching 
fractions. 

For the other background components the authors rely on the MC prediction of the 
differential shapes in and mx- 

The one-dimensional projections of the fitted and measured mx-q^ distribution are de- 
picted in Fig. |3] and the fit has a good statistical significance of 12%. The binning was 
chosen in order to minimize the dependence on the B — Xu I vi signal model. The partial 
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I Kb I X 10^ 



BLNP [4r 



GGOU M 



DGE [35; 



4.37±0.26l°^f 4.41 ±0.26;°^^ 4.46 ± 0.26 



+0.15 
0.16 



Table I The determined values of | V-ah \ of Ref . [53] are listed using the prediction of the partial 
decay rates of Refs. [Ml ESj H] . The uncertainties are experimental and due to theory. 



branching fraction is calculated as 

AB (IpT^I > 1.0 GeV) = X (1 - 6qed) , (5) 

where and correspond to the fitted B — )■ I signal yield and efficiency, and A^tag 
corresponds to the number of tagged B decays. The factor 5qed corresponds the arising QED 
correction which is predicted using the algorithm of Ref. . The systematic uncertainties 
to the partial decay rate were studied by varying the constitution of the B — )• I i^i hybrid 
mix, the shape function moments, and the B ^ X^l I'l contributions. The value of \Vub\ can 
be readily obtained from Eq. ([s]) using the predictions for the partial differential decay rates 
of Refs. [an [351 SI] via 

rub I - Y,^Ac(pr>i.o) w 

where tb denotes the averaged B^ and i?''-meson life time, and the prediction for the 
partial decay rate stripped of the squared CKM matrix element. The determined values for 
I Vub I are given in Table [ij 



B. Study of B ^ Xului from BaBar: Ref. [42] 

The authors of Ref. [12] measure the partial branching fraction of i? — )■ Xu I i^i using 
various kinematic cuts to suppress the predominant B ^ X^lvi background. In particular, 
cuts on the following kinematic variables are studied to separate signal from background 
contributions: the four-momentum transfer of the B meson to the X^ system squared, g^, 
the lepton three-momentum in the i?-meson rest frame, p;, the light-cone momentum of 
the hadronic system obtained the 5-meson rest frame, p"*" = Ex — \'Px\i and the hadronic 
invariant mass, mx- In total 467 x 10^ BB decays measured at the PEP-II B factory 
with the BaBar detector were analyzed using a tagged approach. The recoiling B meson is 
reconstructed using a semi-exclusive algorithm based on hadronic B — )• D'^*'>Y decays, where 
y is a charged system composed of pions and kaons. The kinematic consistency of the B 
meson candidate is checked using the beam-energy substituted mass Eq. ([T]) and the energy 
difference Eq. (|2]). The selection requires a value of /S.E compatible with zero within about 
three standard deviations of the resolution uncertainty. If more than one B meson candidate 
is present in an event, the candidate with the lowest from the summed terms of a 
kinematic vertex fit, the central value of the reconstructed D mass, and the compatibility of 
AE with zero, is chosen. This selection results in at least one B meson candidate in 0.3% 
of all B^B^ and in 0.5% of all B+B' events. 

The B — )■ Xu I vi signal selection requires at least one electron or muon candidate, which is 
identified using a particle identification algorithm, and needs to pass an angular acceptance 
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cut. A lower cut on the three-momentum in the S-meson rest frame, of 1 GeV/c^ 

is imposed in order to suppress background from cascade and r decays. Muon candidates, 
which are consistent with originating from a, J/ip when paired with another charged track 
of opposite charge in the event, are rejected. Further, electrons which are consistent with 
originating from a 7 — > e^e~ conversion when paired with another charged track of opposite 
charge in the event, are rejected. The hadronic system X is reconstructed from charged 
tracks and energy depositions in the calorimeter that are not associated with the tagged B- 
meson candidate. Requiring exactly one electron or muon candidate suppresses B ^ XcIpi 
background which frequently produce a second lepton in cascade decays. In addition, a 
charge correlation cut with the determined charge of the B meson from the tag side, Qb, 
and the charge of the reconstructed lepton, Qi is imposed, such that, Qh Qi < 0. In addition, 
the total charge of the event is required to be zero, i.e. Qtot = Qb + Qi + Qx = 0, where Qx is 
the charge of the hadronic X system. After this selection three main sources of background 
remain: combinatorial background from the tagged side; background from B — )■ X^l^i 
and cascades; background from B — > X„ r ui decays with r — )■ e or /i. The combinatorial 
background is subtracted on the tag side by a maximum likelihood fit of the rriEs distribution. 
The other two sources can be reduced by reconstructing the missing mass from the missing 
four-momentum of the event, i.e. in the T{4S) rest frame 



where p*^ refers to the four-momentum of the tag side, p*x is the four-momentum of the 
hadronic system associated with the B — )■ X„ / z// candidate, and the four-momentum of 
the lepton candidate. For a proper B — > X lui decay the missing mass, m^jgg = p^Ls 
required to peak at zero and it is required that — 0.5GeV < m^j^g < 0.5 GeV. Further, 
a veto similar to Eq. ([s]) is implemented. Finally, if a charged kaon or a Kg is identified 
by a particle identification algorithm in the tracks of the recoiling B meson, the event is 
rejected. The combined selection rejects ^ 90% of the background decays that passed the 
preselection, while retaining about ^ 33% of the signal decays. 

The applied vetos allow the separation of the candidates into two samples: a signal 
enriched sample which includes the candidates that pass all the selection criteria, and a 
signal depleted sample, which contains all candidates that passed the preselection but failed 
at least one veto requirement. The signal depleted sample is rich m B Xclvi events and 
can be used to cross check the background assumptions in the signal enriched sample. Both 
samples are analyzed with an implicit cut in \pi*^\ > 1.0 GeV unless stated otherwise and 
with fits in 

- mx with a cut on mx < 1-55 GeV and mx < 1.7 GeV 

- p'^ and a cut on pj^ < 0.66 GeV 

- two dimensions in mx — (f-, and cuts on mx < 1.7 GeV and > 8 GeV^ 

- two dimensions in mx — q^, and a cut on \pi*^\ > 1-3 GeV in 

- \pi*^\ with a cut on \pi*^\ > 1.3 GeV 

The value of is reconstructed using Eq. ^ and the fit has two free parameters: the 
B — 7- Xu I vi signal yield, and the summed B ^ Xclvi and other background yield. 




(7) 
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Fig. 4 The various projections of the kinematic quantities of the one-dimensional fits to the 
B —7- Xu I vi and background yields are shown. Black depicts the data points, the 
B — >• Xu I vi signal is shown as a histogram shaded in white, B — )• X^ I vi whose actual 
true kinematic values are outside the selected phase-space region is shown in cyan (or 
light grey), and the remaining summed background corresponds to the dark grey 
shaded histogram: (a) mx, (b) pj, (c) , and (d) pi. The bottom row depicts the 
rebinned background subtracted signal yields. 



The B — X^lvi signal is modelled similarly as in Ref. [53], i.e. with a mix of resonant 
and non-resonant decays. The resonant B — )■ nlvi decays are modelled using the model 
of Ref. [19], with the measured model parameter of Ref. [U]. For B rjli^i, B — )■ rj'lui, 
B — >■ plui, and B — >■ uWi the form factor parametrization of Refs. [TU |15] are used. 
The corresponding exclusive branching fractions are fixed at the measured world-averages 
as calculated by Ref. [16]. The non-resonant contributions are modelled using the shape 
function parametrization of Ref. |33j . 

The B — )• Xf-lvi contributions are modelled as the sum of the exclusive B — )■ Dli?i, 
B ^ D* I Ui, B ^ D**li'i, and non-resonant B — )■ D^*^Tili'i decays using the form factor 
parametrizations of Refs. [2S1 ESI 113] . The ratio of the resonant B — )■ D** I vi with respect 
to the summed B — t- D^*"^ I vi plus all other background, i.e. 

\ ^ B{B^D'*lui) /gx 

^^** B(B-!.D* /j/;)+othcr background V°/ 

is determined from a two-dimensional fit in — mx to the signal depleted sample. The 
binning in the signal region was chosen to minimize the dependence on the B — )■ I I'l 
signal model. The results of the one-dimensional and two-dimensional fits for the signal 
and background yields are depicted in Figures |4] and [s] |Vub| can be readily obtained from 
the fitted signal yields and the total inclusive B ^ Xlvi branching fraction using Eq. (|6]). 
The authors use the predicted differential decay rates from Refs. [21 |3l |3ll |35l HI] and Table 
[n] lists the determined values of |Vub| from the two-dimensional fit in and mx with the 
impUcit cut of \pi*^\ > 1.0 GeV. 



C. Summary of inclusive results 



The author of Ref. [16] averaged several recent measurements of \V[ 



ub 



untagged approaches to determine the partial B — X^ I vi decay rate. Table III 



using tagged or 
summarizes 
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Fig. 5 The projections of and mx of the two-dimensional fit in — mx with an imphcit 

cut in \'Pi*^\ > 1.0 GeV are depicted. 





BLNP m] 


GGOU [34] 


DGE [35] ADFR [2 13] 


IKbl X 103 


4.27± 0.23t°:^^ 


4.29 ±0.24^°:^^ 


4.34 ± 0.24 ± 0.15 4.35 ± 0.28^°;^^ 



Table II The result of the determined value of iKibl from the two-dimensional fit of Ref. [12] 
in and mx with an implicit cut in pi > 1.0 GeV are listed. The partial decay rates 
were determined from the predictions of Refs. [21 [3l [Mj [35l I41j . The uncertainties are 

experimental and due to theory. 



the results of the untagged measurements Refs. [3 El EH US] and compares its average with 
the average calculated from the two presented tagged analyses Refs. jlSl |53|. The value 
of iKibl from both approaches are in very good agreement with each other. Further, the 
confidence regions of all theoretical calculations overlap. 

The good agreement of untagged and tagged results indicates that the different experi- 
mental approaches yield consistent results. Further, the good agreement between the var- 
ious QCD based predictions for the inclusive decay rate indicate a sane treatment of the 
uncertainties related to the various approaches and approximations. However, the chosen 
approaches to describe the charmed background and the signal decays lead to some scrutiny: 
Although the used assumptions reflect our best knowledge, it also incorporates some incon- 
sistencies which potentially could have an impact on the determined values of |Vub|- The 
measured inclusive B ^ Xcli'i branching fraction is not identical to the measured exclusive 
contributions from B ^ D I i/i, B ^ D* I ui, and the measured branching fractions of the 
decays into the narrow and broad IP states, B D** I ui, i.e. a gap of about 1.4% branch- 
ing fraction persists. One possible explanation for the presence of this gap would be the 
existence of further strong decay channels of the D** mesons, i.e. the B D** I z/; branch- 
ing fractions is measured only over the reconstruction of the D** via D** D^*h. The 
presence of strong decays with an additional pion or into another light meson would allow 
for further contributions from IP states. Such transitions, however, are not experimentally 
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BLNP 


GGOU 


DGE 


Experiment 


Selection 


ABdO-") 




IVubl (10-^) 




CLEO (1) 


p* > 2.2 GeV 


2.30 ± 0.1 5 ± 0.30 


4.00 ± 0.47 ± 0.34 


3.81 ± 0.45+°-3g 


3.70 ± 0.43+ 


Belle (2) 


p* > 1.9 GeV 


8.47 ± 0.37 ± 1.53 


4.81 ± 0A5*IH 


4.65 ± 0.43+°- 30 


4.66 ± 0.43+ °-25 


BaBar (3) 


p* > 2.0 GeV 


5.72 ± 0.41 ± 0.51 


4.35 ± 0.25!°jJ 


4.17 ± 0.24+°-33 


4.15 ± 0.28+ °-25 


BaBar (4) 


p* > 2.0 GeV 

jmax. > 3 5 2 


4.41 ± 0.42 ± 0.42 


4.48 ± 0.30+°-3' 


— 


4.15 ± 0.28 ± 0.30 


Average (1-4) 


Untagged 




4.41 ± 0.17± 0.32 


4.20 ± 0.19+°-33 


4.16 ± 0.17+°-25 


Belle (6) 


M X X p* > 1.0 GeV 


19.6 ± 1.7 ± 1.6 


4.45 ± 0.27!°;^^ 


4.47 ± 0.27*° ]l 


4.53 ± 0.27 ± 0.15 


BaBar (7) 


M X X q2, pj > 1.0 GeV 


18.0 ± 1.3 ± 1.5 


4.27 ± 0.23!°:^3 


4.29 ± 0.24+° ]{ 


4.34 ± 0.24 ± 0.15 


Average (6) &(7) 


Tagged 


18.7 ± 1.0 ± 1.1 


4.35 ± 0.18+°-2^ 


4.37 ± 0.18+° ]5 


4.42 ± 0.18 ± 0.15 



Table III The determined value of |14ib| from a selected number of untagged and the two 

presented tagged measurements are listed: CLEO (1) corresponds to Ref. [21], Belle 
(2) to Ref. @5], BaBar (3) and (4) to Refs. and |7]. Belle (6) and BaBar (7) are 
the two tagged measurements presented in this talk. The stated uncertainties are 
experimental and from theory. Further, the Table was taken from Ref. |46j . 



verified yet. 

Both analyses adopt a different approach to deal with this issue: Ref. [53] fills the gap with 
broad IP decays and allows the two-dimensional fit to further scale the summed background 
components (consisting of B X^li^i, secondary leptons and misidentified candidates). 
Although a systematic uncertainty is assigned (and the impact on the determined value of 
iKibl seems small), the authors do not specify if the adjustment of the background yields 
is compatible with the measured branching fractions for B — )• Dlui and B — )• D*lvi. 
Ref. [12] fills the gap with broad and narrow IP states. The fit to the signal and background 
yields is further allowed to scale the charmed and other background separately. In addition, 
the authors determine the ratio Eq. (|8]) from their signal depleted sample, what results in 
an overall increase of the IP contributions with respect to the IS* decays and the other 
backgrounds. The obtained value of the partial branching fraction seems fairly insensitive 
to the B Xclvi yields due to the charm suppressing phase-space cuts. 

Further, both analyses make the implicit assumption, that the shape of the other back- 
grounds (e.g. from secondary leptons and misidentified candidates) is correctly described 
in the MC simulation, but their corresponding yields are off. This of course is not very 
satisfactorily. 
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3. EXCLUSIVE DETERMINATION OF |Kb| 

In principle, many exclusive final states could be studied to determine |Vub|, but the the- 
oretical and experimental most promising decay channel is B ^ n I uf. the decay signature 
with a single charged or uncharged pion and a high-energetic lepton allow with the combi- 
nation of the known beam-energy constraints the reconstruction of the missing mass of the 
event, which is associated with the neutrino that eludes detection. This is often combined 
with a tagged approach, i.e. one of the decaying B mesons is reconstructed, what allows 
a more clear separation of background associated with the decay of the second B mesons. 
Such a tagging, however, results in a lower overall reconstruction efficiency, and the three 
analyses presented in this talk use an untagged approach. 

The B ^ irlvi decay rate can be parametrized as a function of a few Lorentz invariant 
amplitudes called form factors. These can reliably be predicted from QCD sum-rules or 
by unquenched lattice QCD calculations and are functions of the four-momentum transfer 
squared of the 5-meson to the pion system, = [ps — Pnf' ■ The value of |Vub| can be 
determined by using the measured partial branching fraction of a given range and the 
predicted partial decay rate. More recent, |Vub| was also determined by a simultaneous fit of 
the unquenched lattice calculations and the measured q^ distributions. In this talk the recent 
measurements of Refs. |37], [29], and |30] are presented. Further, the combined extracted 
value for |V^b| of a fit to the unquenched lattice QCD results of Ref. [12] and the measured 
partial B ^ nlvi branching fractions of Refs. [301 EI] are presented. 

A. Study of S° TV- 1+ u decays from Belle: Ref. [31] 

The authors of Ref. [3Z] study B^ u from 657 x 10^ of BB pairs produced at the 

KEK-II 5-factory and recorded by the Belle detector. The signal decay of B^ — ?■ tt^ ui is 
reconstructed from all oppositely charged leptons and pion candidates, which are required 
to be positively identified by a particle identification algorithm. The lepton is required to 
be either an electron or a muon. Both candidates are fitted to a common vertex and are 
required to possess a significance level of greater than 1%. The missing four-momentum of 
each event is calculated in the T(45') rest frame by 



where ^/s is the average beam-energy, and the sum runs over all charged and neutral particle 
candidates in the event (where {E*,p*) denotes the four-momentum of the z*^ candidate). 
In order to be selected an event is required to have -E^iss > GeV. The missing three- 
momentum of the event is associated with the neutrino and the neutrino four-momentum 
is determined as p* = (IPmisJ )Pmiss) ^o the higher resolution of the three-momentum 
reconstruction in comparison to the energy resolution. To select events compatible with the 
signal decay the total charge of the event is required to be small or equal to three positive 
or negative elementary charges. The cosine of the angle between the combined pion and 
lepton three-momentum and the three-momentum of the initial B meson three-momentum 
in the T(4S') rest frame proofs to be a good discriminator to reject background decays, i.e. 





cos 6 



BY — 



2|Ps| \Py\ 



(10) 
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where py = {Ey,Py) denotes the summed four-momentum of the lepton and pion candidate, 
with mass squared my = py^. The absolute value of the 5- meson momentum is calculated 
from the beam energy as = ^^/s/A — m\. If the lepton and pion candidate are from 
a signal decay, the latter cosine of the angle is required to range from —1 to 1. Final 
state radiation is responsible for a marginal tail into negative cos G^y values, and resolution 
effects shift a small number signal events into the region above 1. Background decays will be 
shifted to large negative values, and uncorrelated lepton and pion pairs will result in a flat 
distribution. In order to suppress background it is required that — 1 < cos G^y < 1. Further 
background separation can be obtained by reconstructing the beam-energy constraint mass 
Eq. ([T]), where the absolute value of the 5-meson three-momentum is calculated as 

\pb\ = \p:+p:+p:\ (h) 

where p* and p* denote the reconstructed and measured three-momenta of the pion and 
lepton candidate boosted into the T(4S') rest frame. The neutrino three-momentum p* is 
inferred from the missing three-momentum of the event. In addition, the energy difference 
Eq. ^ is calculated using 

ei = e:+ e: + e: (12) 

where £"*, and £"* denote the energy of the pion, lepton, and neutrino candidate as 
calculated from the reconstructed three-momenta or the missing three-momentum, respec- 
tively, boosted into the T(45') rest frame. Candidates outside the signal regions, defined 
by \^E\ < 1.0 GeV and rribc/ES > 5.19 GeV, are rejected. To suppress further background 
from continuum decays, cuts on the zeroth, first, and second Fox- Wolfram moments are 
implemented. In addition events with an invariant mass my ranging from 3.07 GeV to 
3.13 GeV are rejected to suppress background from misidentified J/V' — )■ yW"^/^" decays. The 
four-momentum transfer squared from the B meson to the pion is reconstructed using the 
four-momentum of the pion candidate and four-momentum of the 5-meson calculated from 
the beam-energy ^/s and its known mass. The B meson lies in a cone around the Y system, 
and a weighted average over the angular orientation is taken in order to calculate g^. The 
potential signal candidates are binned into 13 bins in q"^, in 16 bins in AE, and in 16 bins 
in rribc/ES- 

The candidates in each bin undergo a further preselection which is tuned to maximize 
the signal to background ratio in each individual bin. The used discriminating observables 
are the angle between the thrust axis of the Y system and the thrust axis of the rest of the 
event; the helicity angle of the lepton-neutrino system; and the missing mass squared of the 
event, i.e. 

"^miss — -^mis ~ IPmissI • (1^) 

The absolute value of the cosine of the angle between the two thrust axes are expected 
to peak near one for signal events, and to be almost uniformly distributed for background 
decays. The cosine of the helicity angle for the lepton-neutrino system, i.e. the cosine of 
the angle between the lepton three-momentum and the three-momentum of the opposite B 
meson in the rest frame of the lepton-neutrino pair, is expected to peak at 1 due to the 
left-handed V — A structure of the weak decay. The missing mass gives an indicator for the 
compatibility of the rest of the event with a decay involving a neutrino. 
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The signal yields are determined by performing a two-dimensional binned maximum 
likelihood fit in AE and rrihc/ES in the plane of the 13 bins, i.e. in total 13 x 16 x 16 
bins are considered. The probability density functions (PDF) describing the signal and the 
various background yields in AE and rribc/ES are obtained using MC simulations. To reduce 
the free parameters, the bins of background yields from h ^ u and b ^ c transitions are 
grouped in a coarser binning, e.g. four and three bins respectively. The distribution from 
continuum background is described by MC, which was reweighed to match the distribution 
of an off-resonance data sample. The continuum normalization is not a free parameter in 
the fit, but fixed at the scaled number of expected off-resonance events. Including the signal 
yields in each bin, there are 20 free parameters in the fit. In particular, 4 yields describe 
the B — )■ Xu I vi background, and 3 yields the background due to B ^ Xclvi decays. The 
resulting number of signal and background yields agree well with the expectations from 
MC simulation studies and the resulting efficiency corrected and unfolded q^ spectrum is 
depicted in Fig. |6| The largest systematic uncertainties of the measurement are due to 
the limited knowledge of the PDF from continuum background, the B^ — > vr" z// signal 
modelling, and the background modelling from B^ — )■ p~ vi decays with p~ — )■ 7r~7r°. 
Other dominant contributions to the systematic uncertainties are due to the performance 
of the particle identification algorithms, the sensitivity on the quality cut on the vertex fit 
probability, the uncertainties due to tracking efficiency, and the composition and modelling 
of the B — )■ Xu I I'l and B ^ Xc.li'i backgrounds. 

The value of |Kib| can be calculated using Eq. ^ and predictions for the (partial) decay 
rates from theory: The unquenched lattice QCD calculations of Refs. [251 115] and the sum 
rule results Ref. [H] predict the value of the form factors at several points at high g^, i.e. 
g^ > 16GeV^. Using a specific parametrization for the form factors, a prediction for the 
partial and total decay rate can be obtained. The values of |Vub| using this approach are 
stated in Table HVl 

This ansatz, however, potentially introduces undesired model dependence due to the 
specific chosen parametrization. This model dependence is somewhat minimized by only 
calculating |Vub| within the g^ range the lattice QCD or sum-rule prediction is considered 
reliable. This, however, is somehow unsatisfying and the authors of Ref. ^30j also explore 
the possibility of a model-independent extraction of iKibl, which makes use of the pro-forma 
model-independent parametrization of Ref. [6] and the predicted normalization from lattice 
QCD of Ref. [19]. In order to perform a combined fit of the parametrization of [6] and 
the points from ^[^, the g^ spectrum, however, has to be transformed into a variable 2, 
which is further discussed in Ref. [18] . The result of this fit is depicted in Figure [7j The 
normalization of the fit can be used to calculate |Vub| and the authors of Ref. [30] quote 

iKbl = (3.43 ± 0.33) X 10"^ , (14) 

where the uncertainty is from both experimental and theoretical sources. 
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Fig. 6 The efficiency corrected and unfolded spectra of — )• vr" vi of Ref |37j is shown. 

The solid curve shows the result of a fit using the parametrization of Ref. [19] . The 
solid histogram shows the prediction using lattice QCD of Ref. pBJ . The dotted 
histograms corresponds to the lattice QCD prediction of Ref. [l9]. The dashed-dotted 
histogram corresponds to the QCD sum-rule prediction of [H]. The dashed histogram 
shows the quark- model prediction of Ref. [21] • All histograms were scaled by | I ^ as 
calculated from Eq. ^ in the cp' range the calculation is deemed reliable. 

gVGeV' iKbl X 10^ 

HPQCD ^ >I6 3.55±0.13l[5j^ 
FNAL [ig > 16 3.78 ±0.141°;^^ 
LCSR [H] < 16 3.64±0.1ll°;^° 

Table IV The extracted values of |Vub| of Ref. [37J obtained using the form factor predictions 
obtained from the lattice QCD points of Refs. |28i 49], and the sum- rule prediction 
of ^4jj. The uncertainties are experimental and from theory. 




Fig. 7 The result of the combined lattice QCD and data points in the parametrization of 
Ref. ^ as presented in Ref. [37j is shown. The variable z is defined in Ref. |18j. The 
determined normalization can be used to calculate |V^b|. 
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B. Study of B ^ Tvlvi and B ^ plui decays from BaBar: Ref. [30] 



The authors of Ref. [3D] study the decays of 5° tt" /+ ui, 5+ 7r°/+ ui, p- 1+ 
and — t- p'^ Z+ z/^ using 377 x 10® BB pairs produced at the PEP-II 5-factory and recorded 
by the BaBar detector. Similar to Ref. f33| an untagged approach is used. The identification 
of the four signal decays require events with at least four charged tracks, the identification 
of an electron or muon using an particle identification algorithm and the presence of one 
or more charged or neutral pions. The lepton candidates are selected from charged tracks 
which have at least a three-momentum of |pf | > 1.0 GeV or \f)l\ > 1.8 GeV for B ^ nlui 
OT B ^ pliJi, respectively, in the rest frame of the T(4S'). To suppress background from 

e~ — e^e~(7), the event as a whole is required to have a minimal longitudinal three- 
momentum such that (z < 0.65 with 

i i 

where the sum runs over all charged particles in the event, and pf and Ei correspond to 
the the longitudinal three-momentum and energy measured in the laboratory frame. All 
tracks used for the reconstruction of the hadron candidate must be rejected by a lepton 
and kaon particle identification algorithm. Neutral pions are reconstructed from pairs of 
photons. The reconstructed vr*^ is rejected if it has not at least a three-momentum of 0.2 GeV. 
Candidates for — )• vr^vr^ and p*^ — )■ 7r"'"7r~ decays are required to have a two-pion mass, 
mTTTT, within the full width of the p-meson mass, i.e. 0.65 GeV < m^^ < 0.85 GeV. To reduce 
combinatorial background, the three-momentum in the centre-of-mass frame of the T(45') 
of one pion candidates has to exceed 0.4 GeV, and further, that the second pion candidate 
has a three-momentum of at least 0.2 GeV. Every hadron candidate is combined with one 
lepton candidates to form a Y candidate, and a vertex fit is performed. A signal candidate 
is required to have at least a vertex fit significance of 0.1%, and a variety of other kinematic 
cuts on the lepton and hadron candidate kinematics are performed to reject background. The 



value of cos 0_By , as defined in Eq. (10), is determined and in order to pass the preselection an 
candidate is required to have —1.2 < cosB^y < 1.1, what is somewhat better than the hard 
cut imposed in the preselection of Ref. [37] , since it reduces the sensitivity due to resolution 
effects. The neutrino four-momentum is reconstructed from the missing energy and three- 
momentum of each event using Eq. ([9]), but in contrast to Ref [37] the missing energy 
and three-momentum are determined in the laboratory frame. Further acceptance cuts are 
imposed to improve the reconstruction of the missing energy. For a correctly reconstructed 
event with a single semileptonic decay, the missing mass squared, m^^gg, is consistent with 
zero, and the failure to detect one or more particles in the event results in a considerable 
tail towards positive values. Candidates are required to have m^-gg/2i?miss < 2.5 GeV, where 
-E-miss denotes the missing energy in the laboratory frame. For the signal extraction the 
beam-energy substituted mass, rribc/ES, and the energy difference, AE, are reconstructed. 
In contrast to Ref. p^, both variables are determined in the laboratory frame, i.e. 



rribc/ES = \l - Pb > (16) 

and 

AE = P^-^e+e_— , (17) 
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where Pe+e- = (-E-e+e-jPe+e-) denotes the measured four momentum of the colhding beam 
particles, and ^/s is the average centre-of-mass energy of the colhding beam particles. The 
S-meson four-momentum is calculated from the measured three-momenta in the laboratory 
frame of the lepton and hadron candidates, and the missing three-momentum of the event 

associated with the neutrino, i.e. pb = Pi +Phadron +pl and ps = {\j rn?B + \'Pb'^ ^Pb)- The 
AE — nibc/ES plane is separated into a fit region, and a sideband region, both within an 
overall signal region, which was chosen to minimize background from B X^lvi decays. 
The four-momentum transfer squared from the 5-meson to the hadron system is calculated 
from the four-momentum of the lepton candidate and the missing three-momentum, as 
determined in the laboratory frame: 

= {Pl + (IPmissI , «Pmiss))^ , (18) 

where a = 1 — AE/E^iss is a scaling parameter which improves the resolution in g^. 

In order to further improve the discrimination between signal and background, a neural- 
network technique based on a multi-layer perceptron is used, cf. Ref. |50]. Seven variables 
are used as input for three neural networks: 1) the angle between the thrust axis of the 
Y candidate and the thrust axis of the rest of the event; 2) the summed absolute values 
of the three-momenta not belonging to the Y candidate weighted with the cosine squared 
of the polar angle of the corresponding track with respect to the thrust axis calculated 
from all tracks not belonging to the Y candidate; 3) cos Q by', 4) "^miss/ (2-E'miss); 5) the 
second normalized Fox- Wolfram moment; 6) the polar angle of the missing momentum in the 
laboratory frame; 7) the helicity angle of the lepton-neutrino pair. The three neural networks 
are trained to either reject continuum, B — )■ Xcli^i, and B — )• Xuli'i background using 
simulated signal and background decays. After the preselection and the filtering through 
the neural networks, the efficiency of selection background from i? — )■ / z/; is < 10^^. The 
efficiency of selecting continuum background is < 10~^, and the efficiency of selecting other 
B — )■ X^lui decays ranges from 0.3 x 10~^ to 0.6 x 10~^. The three neural networks thus 
suppress the background selection of the dominant BB and continuum background by a 
factor of 10^ and 10^. The efficiency of selecting a 5° — t- vr^ vi or i?+ — t- tt^ vi signal 
decay is 1.8 x 10^^ or 1.6 x 10"^, respectively. The efficiency of selection a B^ ui, 
and B~^ — )■ p^l^vi signal decay is 0.3 x 10~^ or 0.8 x 10~^. The rejection of the neural 
networks for simulated continuum events was studied with measured off-resonance data. In 
addition, the understanding of the neural network rejection of charmed semileptonic decays 
was studied with a, B X^l vi enhanced sample, and a 5'^ — ?► D*^ vi control sample. 

At this point, still several signal candidates per event are allowed and the simulation 
reproduces the candidate multiplicity well. In order to reduce the dependence on the sim- 
ulated candidate multiplicity, only the Y candidate with the highest significance from the 
vertex fit is retained. In case for B^ — ^ vr" /+ vi the candidate closest to the central value of 
the 7r° mass is selected. MC studies using these selection rules imply that in about 60% of 
all cases the correct signal candidate is selected. 

The signal yields for i?° vr" /+ z//, _B+ 7r°/+z/;, _B°p"/+i//, and 5+ p°/+z/; are 
determined in an extended maximum-likelihood fit to the three-dimensional AE—mhc/ES~(l'^ 
distributions for each mode, which takes into account the limited knowledge of the signal 
PDFs in the MC simulation. The binning in the AE — mf,c/ES plane was chosen to optimize 
the background and signal shape discrimination by retaining an adequate statistics in all 
bins. Overall 47 bins in AE — rribc/ES are chosen and the range for B — )■ irlui from 
to 26.4 GeV^ is divided into six bins, and the range for B — )■ plui from to 20.3 GeV^ 
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is divided into three bins, resulting in a total number of 282 and 141 bins, respectively. 
In the nominal fit, all four signal yields are simultaneously determined from a fit to the 
measured AE — rrihc/ES ~ distributions from all four channels, taking into account the 
cross feeds, e.g. the signal decay in one data sample may contribute to the background 
in another sample. Furthermore, the isospin relations between the branching fractions of 

—7- TT^ Z+ vi and — t- tt'' z/^, and between B^ p~ z// and B~^ — )■ Z+ ui, are imposed. 
The yields for B — )■ Xu I vi background for i? — t- tt / z// is determined in two regions, i.e. 

< 20GeV, and > 20GeV^. For B p/z/; the B Xuliyi background is fixed, 
due to the lack of discriminative power between the signal and B — ?► Xu I vi decays. The 
background from other BB processes is separated into two yields: the dominant B ^ D* lui 
and 'other' background (consisting of the remaining B ^ X^lvi and other BB processes, 
e.g. background producing secondary or misidentified leptons). In addition, the yields due 
to continuum background is a free parameter. The cross-feeds between the B ^ tt lui and 
the B ^ plui samples is also a free parameter of the fit. In the nominal fit, all background 
parameters which are not fixed, are fitted separately for each signal decay, since the different 
final states lead to different combinatorial backgrounds which are a-priori unrelated. The 
fit significance of 68.5% is excellent and the resulting isospin combined q^ distributions, 
which were unfolded and efficiency corrected, are depicted in Fig. [8j The leading systematic 
uncertainties for B vr / z/^ are due to the uncertainty in the performance of the particle 
identification algorithms which identify the lepton candidates, the limited knowledge on the 
Ki spectrum and reconstruction performance, the tracking efficiency of the detector, the 
photon reconstruction efficiencies, the knowledge of the continuum PDF, and the modelling 
of the B Xclvi and B — )■ Xu I vi background. For B — t- -nlvi the largest uncertainties 
arise from the shape function parametrization and branching fraction used for the modelling 
of the B — )■ Xu I vi background components, the sensitivity to the signal shape obtained 
from the sum-rule prediction of Ref . [15] , and from the continuum background. 

The value of |Vub| can be calculated using Eq. (|6]) and the predictions from lattice QCD 
and sum rules from Refs [HI [151 1211 [32]. They are stated in Table |V| 

The predictions for the differential decay rates of Refs [HI |T21 EH] make use of a specific 
parametrization for the form factors, what introduces an undesired residual model depen- 
dence into the obtained results. The authors of Ref. [30] explore a similar approach than 
Ref. [37] to combine the measured spectra with the unquenched lattice QCD predictions 
of Ref. [281 SS] cind the pro-forma model independent parametrization of Ref. [HI [23] for 
the form factors for a combined fit which determines iKibl- In contrast to the approach 
presented in Ref. [30], the measured spectrum is not transformed into another variable. 
The fit result to the full spectrum and four lattice QCD points of Ref [19] is shown in 
Fig. |9} The fit has a good significance of 35.5% and for the determined value of |Kib| the 
authors of Ref. [SU] quote 

iKbl = (2.95 ± 0.31) X 10-3 , (19) 

where the uncertainty corresponds to the combined experimental and theory uncertainties 
from the lattice QCD points. 
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The unfolded spectra of the isospin combined sample of i? — >■ vr / z?; (left) and 
B ^ plvi (right) of Ref |3U| are shown. Left: The blue dashed line (denoted as BK) 
shows the result of a fit using the parametrization of Ref. |19] . The blue dotted line 

(labeled as BZ) uses the parametrization presented in Ref [H]. The blue 
dashed-dotted line (BCL) uses the parametrization of Ref. [22], and red line (BGL) 
shows the fit result of the parametrization of Refs. [ISJ [23]. The grey shaded region 
corresponds to the uncertainty of the BGL fit. Right: the predictions of Ref. [51] 
(1SGW2), and Ref. [E] (LCSR), where the ISGW2 prediction has been normahzed to 
correspond to the total measured branching fraction, and the LCSR prediction is 
scaled by |Vub| calculated in the region of g2 < 16. 



B^TTlvi gVGeV^ |Kb| x 10^ 

HPQCD [ig >l6 3.21 ± 0.17l°J^ 
LCSR [H < 16 3.63±0.12j:^;^9 

LCSR [33 < 12 3.78±0.13l°;|^ 



B^plvi gVGeV^ I Kb I X 10' 



LCSR [H] < 16 2.75 ±0.24 

Table V The extracted values of | Kib I of Ref. [30] obtained using the form factor predictions 
obtained from the lattice QCD points of Ref. [28] , and the sum-rule predictions of 
[14^ [T5l [32] . The uncertainties are experimental and from theory. 



18 




q2 (GeV^) 

Fig. 9 The result of the combined lattice and data fit of Ref. |30j is shown. The four solid 
magenta points from Ref. [l9] were combined with the full (f' distribution to determine 
I Vub I using the form factor parametrization of Ref. |18| [23] . All lattice points were 
scaled by the determined value of |Vub| • The grey shaded region corresponds to the 

uncertainty of the fit. 



C. Study of S° TT" 1+ VI and B+ -)• t]^'^ 1+ ui decays from BaBar: Ref [29] 

The authors of Ref [23 study the decays of — t- tt^ ui and — t- r/*^'^ ui from 
464 X 10^ BB pairs produced at the PEP-II 5-factory and recorded by the BaBar detector. 
The analysis uses an untagged approach and a reconstructs the B^ tt^ vi and B^ — )■ 
r]^ ^ /"•" ui signal decay by positively identifying a track as a charged lepton by a particle 
identification algorithm. Electrons and muons are required to have a three-momentum 
greater than 0.5 GeV or 1.0 GeV, respectively. Further, for the B^ — )■ vr^ /+ vi decay a positive 
identification from a particle identification algorithm for a pion candidate is required. The r] 
meson is reconstructed using the decays 77 — )■ 77 and rj — )• tt+tt^vt^. The 77' is reconstructed 
by combining 77 candidates with two charged pion candidates, i.e. through the decay chain 
of 7]' — )■ Tj-K^-K^ . Each possible combination of a lepton and a hadron candidate form a signal 
candidate, referred to as a y candidate in the following. For each Y candidate a vertex fit is 
performed and candidates are required to have a fit significance of larger than 1%. Further, 
angular acceptance cuts on the reconstructed lepton and hadron candidates are imposed 
to make sure that the events were reconstructed entirely within the detector. To avoid 
background from J/ip decays, Y candidates are required to be outside the J/psi 

mass window. In addition, further kinematic cuts on the lepton and hadron four-momenta 
are imposed to improve the signal to background separation. The missing four momentum 
of the event is calculated using Eq. (|9]). The four-momentum of the neutrino is inferred from 
the missing three-momentum of each event, i.e. = (IPmissI ' Pmiss)- The four-momentum 
transfer squared from the B meson to the hadron system is calculated as 

= {PB -Phadron)^ , (20) 

where Phadron denotes the four-momentum of the reconstructed hadron candidate in the 
laboratory frame. The B meson lies in a cone around the Y system, and a weighted average 
over the angular orientation is taken in order to calculate q^. The spectra is further 
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B^TTlDi gVGeV^ |Kb| x 10^ 

HPQCD 28 >T6 3.28±0.20i;]j^ 
FNAL [49] > 16 3.14 ± 0.18l°;^^ 
LCSR [32] < 12 3.70±0.1ll°j^ 

Table VI The extracted values of |Vub| obtained using the form factor predictions obtained 
from the lattice QCD points of Refs. [MlllH], and the sum-rule prediction of [32] . 
The uncertainties are experimental and from theory. 



unfolded to correct for resolution effects. In contrast to Ref. [30] a cut-based approach 
similar to Ref. [37] is used to boost the signal to background ratio for each bin. The 
following cuts are optimized: the angle between the thrust axis of the Y candidate and the 
thrust axis of the rest of the event is optimized for each region; the polar angle of the 
missing momentum in the laboratory frame; the missing mass squared divided by the missing 
energy, m^jgg/ (2i?niiss); and the helicity angle of the lepton-neutrino pair. On average, about 
1.14 candidates are observed in each selected event, and for events with multiple candidates 
only the candidate with the larges value of the lepton-neutrino helicity angle is kept. The 
signal selection efficiency varies between about 8% to 15% for — ^ ui (depending 

on the range), and between about 1.5% and 2.6% for — >■ rjl^ui decays in case the 
Tj was reconstructed via two photons. In case the rj was reconstructed via three pions, the 
reconstruction efficiency drops to about 0.6%. The reconstruction efficiency of B~^ — )■ i]' /+ ui 
is also about 0.6%. The 5° — )■ vr^ ui signal selection efficiencies are considerably higher 
when compared to the selection efficiencies of Ref. [30], but so are the background selection 
efficiencies. 

The signal yields are determined by a two-dimensional binned maximum likelihood fit in 
AE and mt,c/ES, as defined in Eqs. (17) and (16). Only candidates with |A£'| < 1.0 GeV and 
TTT'bc/ES > 5.19 GeV are retained. The 5° — )■ ui background is grouped into five cate- 
gories: 1) B ^ Xu I vi originating from the same B as the signal candidate; 2) B ^ I vi 
originating from the opposite i?; 3) BB background originating from the same B; 4) BB 
background originating from the other 5; 5) continuum background. The B^ — )■ tt^ Z+ vi 
distribution is divided into 12 bins for the signal yields, and two bins for each of the five 
background yields. The bin widths of the background components is chosen in such a way 
that both bins retain adequate statistics. The PDFs for the AE and mi,c/ES distributions 
for signal and background are determined from MC simulations. All 12 signal and two times 
five background yields are determined simultaneously and the resulting unfolded and effi- 
ciency corrected spectrum is depicted in Fig. 10 For B^ rjl^ vi the B — )■ X^ I vi and 
BB background categories from the opposite and same B meson as the signal candidate 
are merged into a single yields. The B^ rjl^ vi distribution is divided into three bins 
for the signal yields, and one bin for the background components. The 3 signal, the BB 
background, and the continuum yield are simultaneously determined in the fit and the back- 
ground contributions from B — )• X^ I i^i are fixed to the MC prediction. The distribution 
for B^ — )■ ?7 /"•" z/; is as well depicted in Fig. 10 The leading systematic uncertainties in both 
measurements are due the modelling of the B — )■ X^ I i>i background and detector effects. 

The value of | Kibl can be calculated using Eq. ^ and the predictions for the partial decay 
rates from lattice QCD and sum rules from Refs [2HI |32l |19]. The corresponding values for 



iVubl are stated in Table VI 
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Fig. 10 The unfolded spectra of S° vr" /+ ui (left) and ^ r]l+ ui (right) of Ref [29] 
are shown. Left: the solid black and green curves show the result of a fit using the 
form factor parametrizations of Refs.[18i 123] and [19]. The red wide and narrow 
dashed lines corresponds to the prediction of Refs. [28] and [l9|, and the blue dotted 

line corresponds to the sum-rule prediction of Ref. [32]. The lattice QCD and 
sum-rule predictions were scaled with the corresponding values of |Vub| squared, cf. 
Table [VTl Right: the blue dotted line corresponds to the sum- rule prediction of 
Ref. [13], normalized to the measured branching fraction. 



D. Summary of exclusive results 



The author of Ref. [46j also averaged several recent measurements of the partial and 
total B ^ Tclui branching fraction using either untagged or tagged methods. Table VII lists 
the obtained averaged, which are in very good agreement in the high region, indicating 
that both experimental methods yields consistent results. The agreement in the low 
region and the total branching fraction are less satisfactorily, hinting a systematic difference 
between both approaches. The background in the low region is dominated hj B ^ Xcl i^i 
contributions, and one possible explanation for the observed difference might be due to 
different model assumptions of the charmed background. 

Table VIII lists the obtained values of | Vub| calculated from the averaged partial branching 
fractions using the predictions of Refs. [HJ [28] . The large uncertainties on the sum rule and 
lattice QCD prediction for the differential decay rate (stripped from the CKM matrix element 
squared) results in a compatible result for both g^ regions. 

Table IX compares the value for \V^ub\ from the untagged averages of the low and high 
g^ regions, as calculated by using the predictions of Ref. [H] and Ref. [2H], with the result 
of the combined data and theory fits of Refs. [37] and [30], which use the parametrization 
of Refs. [iHl |23|. Needless to say, such fits fill an important gap and their implications 
are interesting. The determination of \Vah\ from decay rates obtained by fits of model 
parametrization to sum rule and lattice QCD predictions neglects available information in 
the g^ spectrum. Using a pro-forma model independent parametrization of the form factors 
allow a combined evaluation of the lattice QCD points at high g^ region, and the entire 
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Btot 


AB(q2 < 12GeV ^) 


AB(q2 > 16GeV ^) 


Average untagged 


1.44 ± 0.03 ± 0.04 


0.84 ± 0.02 ± 0.03 


0.36 ± 0.03 ± 0.02 


Average tagged 


1.32 ± 0.08 ± 0.03 


0.69 ± 0.05 ± 0.02 


0.37 ± 0.04 ± 0.01 



Table VII Summary of tagged and untagged averages for the measured B" — t- vr^ vi 

branching fractions, as determined by Ref. [l6]: For the tagged average the results 

of Refs [U [21 Uni SO] were used. The untagged average was calculated from the 
measured total and partial branching fractions of Refs. |29 | 1301 [37j . i.e. the three 

presented measurements. 



Theory 


Experiment 


range (GeV^) 


AB= (10-^*) 


Af''(ps-i) 


IVublMIO-^) , 


LCSR 


untagged average 


0-12 


0.84 ± 0.04 


4QQ+1.01 

^■^^-0.95 


3.69 ± 0.08+°-^^ 


tagged average 


0-12 


0.69 ± 0.05 


4QQ+1.01 

^■""-0.95 


3.37 ±0.13!°:" 


HPQCD 


untagged average 


16-26.4 


0.36 ± 0.04 


2.02 ± 0.55 


3.41 ± 0.10!°;^^ 


tagged average 


16-26.4 


0.37 ± 0.04 


2.02 ± 0.55 


3.47 ± 0.19!°-^° 



Table VIII Summary of the tagged and untagged values of | V^b I determined using the 

predictions of Ref. [32] (LCSR), and Ref. [28] (HPQCD), as determined by 
Ref. [l6]: The untagged averages correspond to the presented results 
Refs. [2S1 [301 [37] and the tagged average are calculated from Refs [D 13 dHl HO] . 
The obtained value of |Vub| from the presented measurements are highlighted in 

magenta. 



measured spectrum. Furthermore, combining experimental and theoretical information 
to to determine the non-perturbative shape parameters of the form factors of Refs. [IHl [23] 
minimizes the uncertainty on the predicted differential decay rate, and consequently on \ Vnh\- 
This is true under the premise that neither theory (e.g. through underestimating theoretical 
uncertainties), nor experiment (through e.g. underestimation of systematic effects or a wrong 
prediction for resolution effects) bias the outcome of the combined fit. These of course are 
non-trivial assumptions for both and lead to some scepticism towards the obtained value 
of |Vub|- However, if this premise really is violated one way or the other, the outcome of 
the determination of \V^\^\ based on the fits of model functions to sum rule or lattice QCD 
predictions is biased as well. 



Fig. 11, finally, shows the fit of the parametrization of Refs. [TSl [23] to the measured 
spectra of Refs. |37J and [3D] and the four lattice QCD points of Ref. [IH] (shown in bold 
magenta markers). The fit assumes that the systematic uncertainties of both spectra are 



100% correlated and was performed by Ref. [3l]. The obtained value of |Vub| is given by 

iKbl = (3.25 ± 0.12 ± 0.28) x 10^^ , (21) 
where the uncertainties are experimental and due to theory. 
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Theory 


Experiment 


range (GeV ^) 


AS" (10-") 


Ai:''(ps-') 


1 |v„b|Mio-^ 


LCSR 


untagged average 


0-12 


0.84 ± 0.04 


4 00+'-'" 


3.69 ± 0.08+°-^^ 


HPQCD 


untagged average 


16-26.4 


0.36 ± 0.02 


2.02 ± 0.55 


3.41 ± 0.10+2^^ 


FMAI -M\\ C 


BaBar PRD:83:032007 


0-26.4 


1.41 ± 0.09 


— 


2.95 ± 0.31 


riNrtL IVilLL. 


Belle PRD:83:071101 


0-26.4 


1.49 ± 0.08 


— 


3.43 ± 0.33 



Table IX The determined value of | V-^h \ calculated using the predicted decay rates at low and 
high q^, i.e. Ref. [3^ (LCSR) and Ref. [2H] (HPQCD), are compared with the result 
of the combined fit of the whole measured spectrum of Refs. [3^ and [30], and 
the lattice QCD predictions of Ref. jl9j (highlighted in blue). 




0-26.4 3.25 ±0.12 ± 0.28 



Fig. 11 The fit of Ref. [31] combining the experimental input of Refs. [301 137j with the lattice 
QCD calculation of Ref. [l9] using the parametrization of Refs. [HI |23j is shown. 
The lattice QCD points used in the fit are shown as filled magenta triangles and are 
scaled by the obtained value of ll^ibl squared. 



4. CONCLUSIONS 



Comparing the averaged tagged inclusive result of |Vab| in Table III, e.g. using the BLNP 
prediction, 

iKbl = (4.35 ± 0.19^^1^) X 10"^ , (22) 

with the combined data and lattice QCD fit results using Refs. [SHI EZ] and Ref. [IH], i.e. 

iKbl = (3.25 ± 0.12 ± 0.28) x lO'^ , (23) 

illustrate the tension between the values of |Kib| obtained from inclusive and exclusive ap- 
proaches: Both values differ by about ^ 2.5 cr, if considered being completely uncorrelated. 
The source of this tension is a long standing issue and interesting due to the fact, that the 
theoretical and experimental methods of both approaches for the determination of the par- 
tial branching fraction and |Kib| are often considered independent. Recently, the author of 
Ref. [27] (see also Ref. [21]) proposed that the presence of right-handed currents could bring 



the results obtained from inclusive and exclusive results into perfect agreement. Fig. 12 



depict the determined values |Vub| quoted in Eqs. (22) and (23) as a function of the right- 



handed admixture en, as proposed in Ref. [27]: Both results overlap with an admixture of 
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-0.4 -0.3 -0.2 -0.1 0.1 0.2 0.3 

Fig. 12 The averaged value of iVubl from the BLNP result (green) and the combined data 
and theory fit to Refs. [301 EZ] R^f. [l9] (yellow) are shown as a function of the 

right-handed admixture e/j. 



about en ^ —0.20. A direct search for right-handed contributions can be realized in decays 
with a non-trivial V — A and V + A admixture (e.g. B — )■ plvi). 

Another potential explanations for the tension between the inclusive and exclusive val- 
ues of iKibl niight lie in the model assumptions for the shape function. Their uncertainty 
estimates raise some flags (see e.g. Ref. [52j) and a model independent determination of 
the shape function would be very desirable. The authors of Ref. [20] reported some recent 
progress by determining the shape function with absorbed l/m^ corrections from B ^ Xg'^ 
decays using the proposed methods of Ref [B]. Fig. 13 shows the determined functional 
form. 

The situation of the various different approaches of modelling the B ^ Xclvi background 
is also somewhat unsatisfying. Even after a decade of having two B factories running, we 
still don't know what exactly consist more than ^ 10% of the total B Xf-lvi branching 
fraction. Although many analyses reject or suppress these background components consider- 
ably, this lack of knowledge limits our understanding of inclusive and exclusive B — )■ X„ / z/; 
decays. The precise determination of the charmed decays should be a priority to everyone 
who wants to measure charmless decays below the endpoint of the lepton spectrum, in the 
low region, or at high mx- At the time being, it is this poor understanding of higher 
charmed resonances that pushes inclusive analyses into regions of phase-space where shape 
function effects play an important role and maybe the discrepancy between the inclusive 
and exclusive values of |Vub| does not hint towards new physics, but just poorly understood 
QCD effects. 

With potentially two super _B-factories being built, we hopefully will be able to investigate 
and eventually resolve this conundrum. 
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k [GeV] 

Fig. 13 The B — )• Xgj shape function from Ref. [20] is shown. It was obtained by a fit to 
three measured B — t- Xg 7 photon energy spectra and a pro-forma model-independent 
expansion of the shape function into a set of orthogonal basis functions. The 
non-perturbative physics of the shape function is described by the expansion 
coefficients and the different coloured curves and confidence regions correspond to 
the fit with two (cqi), three (C012), four (C0123), and five (C01234) basis functions. The 
confidence region cover the experimental uncertainties only. 
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